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Spectral properties of oscillatory and non-oscillatory a 2 -dynamos. 
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The eigenvalues and eigenfunctions of a linear o 2 -dynamo have been computed for different spatial distributions of an isotropic cv-effect. 
Oscillatory solutions are obtained when a exhibits a sign change in the radial direction. The time-dependent solutions arise at so called 
exceptional points where two stationary modes merge and continue as an oscillatory eigcnfunction with conjugate complex eigenvalues. 
The close proximity of oscillatory and non-oscillatory solutions may serve as the basic ingredient for reversal models that describe abrupt 
polarity switches of a dipole induced by noise. 

Whereas the presence of an inner core with different magnetic diffusivity has remarkable little impact on the character of the dominating 
dynamo cigenmodes, the introduction of equatorial symmetry breaking considerably changes the geometric character of the solutions. 
Around the dynamo threshold the leading modes correspond to hemispherical dynamos even when the symmetry breaking is small. This 
behavior can be explained by the approximate dipole-quadrupole degeneration for the unperturbed problem. 

More complicated scenarios may occur in case of more realistic anisotropics of a- and /3-effect or through non-linearities caused by 
the back-reaction of the magnetic field (magnetic quenching). 
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1 Introduction 

Oscillatory or reversing magnetic fields driven by a flow of an electrically conducting fluid are a well 
known astrophysical phenomenon. Thus the solar magnetic field regularly oscillates on a characteristic 22 
yrs time scale whereas the dominating dipole component of the Earth's magnetic field irregularly changes 
its orientation every few 100 kyrs, conducting a so called reversal. Sign changes of the magnetic field 
also have been observed in the Cadarache von-Karman-Sodium (VKS) dynamo, and depending on the 
difference of the rotation rates of the two flow driving impellers various regimes with oscillatory and/or 
chaotic behavior can be obtained (jBerhanu et all 120071 ). Whereas the time scale of the solar cycle can be 
reproduced using simple afi-dynamo models, an explanation of the observed equatorial migration requires 
further assumptions, e.g. meridional flow or negative radial shear at the surface. Regarding the geodynamo, 
three dimensional simulations of the magnetohydrodynamic equations as well as mean field models have 
been able to reproduce essential features of the Earth's magnetic field (jHovng et all [2002J, IChristensen 
201 ll ) , but others, for example the very nature of the reversal mechanism or the large variation of reversal 
rates (jMcFadden and Merrill Il986l ) still are unsolved issues. A prominent feature of this paleomagnetic 
reversal frequency distribution is the occurrenc e of a few very long periods 20Myrs) during which no 
reversal occurred at all (so called superchrons; Harland et al . 1982). Another surprising but less known 
property are deviations of the distribution of int er-reversal time periods from an ideal Poisson distribution 
(jCarbone eiail I2006L [1 orriso- Valvo et all l2007n indicating long term correlations in the reversal trigger 
mechanism(s). 

In order to reliably disentangle internal reversal trigger mechanisms (e.g. intrinsic changes of the fluid 
flow pattern in the liquid outer core) from external sources (e.g. precession or changes in the heat flux 
through the core mantle boundary) numerical simulations are required that cover sufficient long time 
periods to allow for a large number of reversal events. In addition, numerical models should also be 
capable to incorporate the mechanisms that trigger individual reversals and reproduce general reversal 
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characteristics such as duration and field geometry during the actual reversal. Mean field models are 
ideally suited to execute long term simulations of the (mean field) induction equation because they are 
computational cheap. The price for this capability is a strong simp lification by parameterizi ng induction 
effects of the turbulent small scale flow e.g. in terms of the a-effect ( Krause and Radler 1980l ). Hence, any 
internal changes of the small scale turbulence (e.g. changes in statistical properties of the convection caused 
by gradual growth of the solid inner core) are suspected to encroach/expand into the comparably simple 
mean field coefficients. Nevertheless, mean field modelling of reversing magnet ic fields has been remarkably 
successful ( Hoyng et al. 20021 . Stefani and GerbetdliooalFischer 6^0!! 12009). The present work is based 



on th e idea that a polarity change is part of an oscillation of the dominant dipole mode (jStefani et al. 



20071 ). A consequential model for irregularly occurring reversals requires a proximity of an oscillatory and 
a non-oscillatory branch so that a single transition between both states might be induced by noise. This is 
substantiated by the fact that polarity reversals in numerical dynamos are generally found in intermediate 
parameter regions between stable (dipolar) dyn amos with small fluctuations and highly fluctuating and 
unstable dynamos ( Olson and Christensen 20061 ) . 

The mean flow in the Earth's fluid core most probably is weak so that usual geodynamo models are 
based on an a 2 mechanism where an a-effect (re-)generates the poloidal field from the poloidal field 
and vice versa. A general requir ement for the occurrence of oscillating eigenmodes in a 2 -dynamos is 
Va 7^ ( Radler and Brauer 19871 ). However, for a long time it was assumed that the leading eigenmode 
in simple a 2 -dynamos is non-osci llatory and dominan t oscillatory solutions are a curiosity that requires 
a rather elaborate configuration (jRiidiger et aZ.I I2Q03I ) . Recently, it has been discovered that oscillatory 



a>- 



dynamos are quite common when t he radial profile of a spherically symmetric, isotropic a exhibits a 



sign change ( Stefani and Gerbeth 20031 ) and a stringent math ematical treatment o f this model yields very 
general conditions for the occurrence of oscillating solutions ( Giinther et al. 2010l ). Meanwhile oscillating 
solutions have also be en found in direct numerical simulations of thermal convection in a spherical she ll 
( Schrinner et al. 2011 ) or in a spherical wedge geometry with random helical forcing ( Mitra et al. 2010l ). 
S uch oscillatin g a 2 -d ynamos might provide an alternative approach for the solar dynamo, e.g. the model 
of iMitra et all (j201(]h exhibits equatorward migration without the requirement of meridional circulation 
or negative radial shear (which in the sun only occurs in a rather thin layer located close to the surface). 

Here, we examine the behavior of axisymmetric eigenmodes generated in a two-dimensional kinematic a 2 - 
dynamo with isotropic a-effect. The focus of our examinations is on the spectrum of eigenvalues in terms of 
growth rates and osc illation frequen c ies for two different radial profiles of the a-effect. For this purpose we 



pick up the model of lGiesecke et al\ (j2005al ) where a sinusoidal radial a distribution is assumed. We show 
that for this a-distribution the leading modes around the dynamo threshold are oscillating independent of 
inner core condu ctivity or latitudinal p erturbations that would provide an equatorial symmetry breaking. 



In extension to iGiesecke et al\ (j2005al ). where only the two dominating oscillating modes around the 
onset of dynamo action were discussed, here we investigate the emergence/disappearance of oscillating 
modes at so called exceptional points in more detail and briefly identify the fundamental differences in the 
radial structure of oscillatory and non-oscillatory eigenmodes. The work aims at building a bridge between 
original mean field models that have been developed in the sixties of the last Century, and recent results 
from observational and numerical approaches demonstrating the valuable impact of mean field theory for 
the development of a geodynamo reversal model. 



2 Equations and Method 

The basic assumption of the mean field approach is a separation of the velocity field u and the magnetic 
field B which are split into mean parts (u), {B) and fluctuating parts u' , B' according to 



B = (B) + B' and u = (u) + u' 



:1) 



where (•) represents an appropriate space- and time average so that the Reynolds averaging rules apply. 
Then the temporal development of the mean magnetic field {B) is described by the mean field induction 
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equation 



d t (B) = V X ((«) x (B) + S - 77V x (B)) 



(2) 



with the turbulent electromotive force (EMF) E = (u' x B') that describes the average induction action 
of (unresolved) small scale flow perturbations. S results from the interaction of the fluctuating velocity 
field ul with the mean magnetic fiel d (B) and is formally rep resented by a linear functional of (u), (B) 
and the statistical properties of u' ( Krause and Radlerl 1980) • Assuming that the mean quantities vary 
only slightly around a certain space-time point, only contribution from a certain neighborhood must be 
taken into account and the EMF can be written as a Taylor expansion: 



V x B'\ 



aij (Bj) + fcjk 



d(B k ) 
dxj 



+ ... 



(3) 



where the coefficients and fiijk depend on the properties of the turbulent fluctuations u'. In general, the 
turbulence is anisotropic so that a- and /3-effect are described by tensors of 2nd and 3rd rank, respectively. 
Here we restrict ourselves to the simplest case of isotropic turbulence so that a and j3 are given by 



a 



'j 



aoSij, 



(4) 
(5) 



so that a and j3 are prescribed by scalar quantities. Note that for a non-vanishing a-effect the turbulence 
additionally must be non-mirrorsymmetric. With the definition of the turbulent diffusivity 



7?T = V + & 



(6) 



and vanishing mean flow {u) = the mean field induction equation simplifies to an equation describing 
the temporal development of an a 2 -dynamo: 



d t B = V x (aB - r/T V x B) 



(7) 



(the brackets around the mean field B from now on are dropped for simplicity). Equation ([7]) is linear in 
B and the ansatz B = B(r)e yt leads to a linear eigenvalue problem 



MB = jB 



(8) 



with the matrix A4 containing the dynamo operator from the right hand side of ([7]), and the eigenvalue 
7 = o+iv (gro wth rate a and frequency v ). For uniform distributions of a and ryx (semi-) analytic solutions 
are known (e.g. Krause and Radlerl . 1980l ). For more complex flows and spatial distributions of a and n a 
nume rical solution of the eigenvalue problem is required (e.g. Roberts and Stixl 1972 . Dudley and Jamesl 
1989). He re we solve (El) i n a sp here surrounded by a non-conducting vacuum applying the method pre- 
sented in Schrinner et al. (l20inh . The approach is based on the biorthogonality of the electric current 
j = /^q X V x B and the vector potential A and explicitly utilizes an expansion of the magnetic field 
B into (analytical known) free decay modes. The method allows a very fast computation of eigenvalues 
and eigenfunctions for nearly arbitrary spatial distributions of a and/or ryx- Further we restrict our ex- 
amination to the axisymmetric field which is dominati ng for isotrop ic a and /3. Higher azimuthal modes 
will be important in case of anisotropic a co efficient (iTilgnerl 120041) . although in case of anisotropic t/t 
axisymmetric modes become again dominant (jElstner and Riidiger 2007|). 
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eigenmodes for (X o =0 



eigenmodes for a =15 




Figure 1. (colour online) Central panel: growth rates for an a 2 dynamo with a = ao cos Note the transition between different types of 
degeneration when going from ao = (free decay) to ao 2> c*crit- The small panels surrounding the central figure present the geometric 
structure of the eigenmodes for ao = (left hand side) and for ao = 15 (right hand side). The colored contours denote the toroidal part 
and the arrows denote the poloidal part. Note the separation of toroidal and poloidal components in the free decay case (ao = 0, left). 

3 Results 

We start with a full sphere (embedded in vacuum) where a is created from turbulent motions in a convective 
layer subject to rotation. Then the EMF contains a term parallel to the mean magnetic field £ oc (g-fl)B 
where g points in the radial direction and fi is parallel to the rotation axis (z-axis) so that (g-Cl) oc cos 1 !?. 
Thus, the a-effect is maximum at the poles and vanishes at the equator: 



a = ao cos #. 



(9) 



This a model has been examined e.g. in iRobertsI (Il972h . Here we revive these results and additionally 



present eigenvalues and eigenfunctions of higher order modes. The resulting spectrum of the leading dy- 
namo eigenmodes in terms of the growth rates versus the amplitude ao is shown in figure Q3 

All eigenvalues are real, i.e. all eigenfunctions are non-oscillatory modes. Due to the strict antisymmetry 
of a with respect to the equator the parity of the eigenmodes remains preserved for increasing a. That is, the 
a-term in the induction equation couples toroidal and poloidal modes in a way that the dynamo eigenmode 
remains antisymmetric (in the following we call this " dipolar- like" ) or symmetric (in the following we 
call this " quadrupolar-like" ) with respect to the equator. At ao = (free decay) an exact degeneration 
between "consecutive" modes is obtained (as indicated on the left side of figure [I]). The degenerated 
eigenfunctions are always a pair of a purely toroidal and a purely poloidal mode which differ in multipolar 
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degree by Al = 1. The degeneration vanishes for ao > and around the dynamo threshold a slight 
predominance of the dipole mode occurs. The critical a for the onset of dynamo action of the dipole 



mode (a.Q = 7.645) is close to the v a lue fo r the quadrupolar mode («q = 7.813). Both values agree 
with the results obtained by Roberts! ( 19721 ) (ag rit = 7.637 for the dipole mode and ag nt = 7.803 for 
the quadrupole mode). For a a cr it all eigenfunctions again approach a twofold degenerated state that 
consists of a pair of degenerated modes with similar geometric structure, but opposite equatorial symmetry 
(see right side in figure [[]). The approximate degeneration between dipolar-like and quadrupolar- like modes 
is the result of a symmetry in the equations describing the a 2 -dynamo which is exact in case of perfect 
conducting boundary conditions. In that case dipolar eigenfunctions an d quadrupolar eigenfunctions are 
adjoints of each other and thus have the same (conjugate) eigenspectrum ( Proctor 19771 ). The degeneration 
is approximat ely retained in case of in sulating boundary conditions so that the eigenvalues remain close 
to each other ( Dobler and Radler 19981 ). 



R adial dependen c e of a Next, we introduce a radial dependence for a taken from the model presented 



m 



Giesecke et al. ( 2005a! ) . This model utilizes an idealized parameterization in the radial direction that 



is based on local simulations of rotating magnetoconvection assumi ng conditions roughly s uitable for the 
geodynamo (weak stratification, fast rotation, strong magnetic field; Giesecke et al. , 2005b! ) . Qualitatively 
these simulations show that on the northern hemisphere a is negative in the lower part of the liquid outer 
core and positive in the upper part, whereby the zero occurs roughly in the middle of the convective 
instable layer. The idealized radial a-profile that incorporates these properties is given by 



a(r) = aocosi?sin ( 2tt 



Rh 



(R 



out 



Rh 



(10) 



where the outer radius of the spherical domain is fixed to i? ut = 1 an d Rm represents the radius of an 
inner co re (see figure [2]). A change in the radial profile was also considered in the first planetary mean field 
models (jSteenbeck and Krausdll969h . However in that work a was related to the radial derivative of the 
turbulence intensity u^ ms and the authors exclusively looked for non-oscillating solutions. 





Radius 



Radius 



Figure 2. (colour online) Radial and latitudinal distribution of the a-effect as given by jlOt . Left: i?; n 
sponding to the actual size of the Earth's inner core. 



0; Right: i?; n = 0.35 corre- 



Here, we start assuming no inner core, i.e. Ri n = 0. The spectrum for the leading axisymmetric modes is 
shown in figure El In contrast to the simple a-profile examined in the previous section, oscillatory solutions 
(indicated by dotted curves) are obtained for this new radial a profile. In particular, around the dynamo 
threshold the dominating eigenmodes are oscillatory dipolar-like (blue curve) and quadrupolar-like (red 
curve) solutions. The oscillating solutions appear/disappear at so called exceptional points (EP) where 
two stationary modes coalesce and continue as two oscillatory modes with conjugate complex eigenvalues, 
i.e. with the same growth rate and positive and negative frequency, respectively (the right panel of figure [3] 
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shows only the positive frequency branch) . Regarding the appearance of oscillatory modes at a 



= (dotted 




curves on the left panel of figure [3]) it is necessary to stress that exactly at a = (i.e. for free decay) these 
modes are stationary and degenerated, i.e. a = is an exceptional point for these modes. The spectrum 
is far more complex than in models without any radial structure of the a-effect and which is apparent 
by the confusing structure of the eigenvalues with a variety of level crossings and/or exceptional points 
where oscillatory solutions constitute or vanish. However, for a 3> a cr ;t predominantly stationary modes 
are observed. 



Decomposition in free decay modes The change of the temporal character of the dynamo eigenmodes 
from non-oscillatory to oscillatory behavior is attended by a change in the radial structure. This can be 
seen by means of a decomposition of a dynamo eigenmode in terms of the leading free decay modes (which 
also represent the basis utilized in the numerical scheme). The dominating contributions in dependence 
on the magnitude of the a-effect are shown in figure HI For both symmetry classes the oscillatory regime 




Figure 4. Contribution of characteristic free decay modes to the first dipolar-like dynamo eigenmode (left panel, eigenmode 1 in figure[3} 
and the first quadrupolar-like dynamo eigenmode (right panel, eigenmode 2 & 3 in figure [3j) . 



is characterized by domination of modes with radial wavenumber n = 1 whereas stationary solutions are 



i, 2012 3:2 Geophysical and Astrophysical Fluid Dynamics giesecke'osc'alpha2 

7 

characterized by a domination of (n = 0) modest- Another feature shown in the decomposition of the 
dynamo eigenmode is the increment of the multipolar degree of the dominant contribution with increasing 
a. 



Inner core The influence of an inne r core with f inite electrical conductivity had been examined by 
Hollerbach and Joned (|l99l Il995h and lWichtl (|2002h with surprisingly opposite conclusions. Whereas the 
latter only found little influence o f finite inner core conductivity on temporal behavior of the dipole field, 
Hollerbach and Joned ( 19931 . 19951 ) concluded from simulations that the finite conductivity of an inner core 
has a stabilizing impact and reversals could only occur if the field in the fluid outer core exhibits a large 
and long lasting fluctuation that allows the field to reverse throughout the inner core as well. 

Figure [5{a) shows that the spectrum of the eigenmodes is only slightly changed when an inner core is 
considered with a radius Ri n = 0.35 and a uniform diffusivity within inner and outer core (t/t = !)• The 
a-effect is prescribed by the distribution (|10p with a = for r < R[ n . The results for the same size of the 
inner core, but with a reduced diffusivity for r < R[ n , is shown in figures [5^b) and (c). The introduction of 
an inner core without changing the diffusivity distribution has remarkable little influence on the spectrum 
(figure [5k,) whereas a reduced inner core diffusivity for sufficiently small a results in a split-up of the 
oscillatory modes into two stationary branches (figure [5b and c). However, around the onset of dynamo 






Figure 5. Growth rates versus oo for a modified o-profile given by ilOh An inner core is considered with an (earth-like) radius R ln = 0.35 
with a magnetic diffusivity r] COIC = 1,0.5,0.1 (from left to right). Solid black curves represent non-oscillatory solutions and grey/dotted 
curves show oscillatory solutions. 



action the typical pattern of the leading dynamo eigenmodes is hardly changed with increasing inner core 
conductivity. For example, the reduction from ^core = 

1 to 

??core — 0.1 leads only to a small decrease of the 
critical a-magnitude from a cr ;t ~ 16.9 to a cr it ~ 16.5. Furthermore, the temporal behavior of the leading 
modes is not changed within this regime, i.e. the leading modes remain oscillating around the onset of 
dynamo action independent of the inner core conductivity. In particular for a ^> a cr it no changes can be 
observed in the behavior of the leading modes. The main impact of an enhanced inner core conductivity 
is manifested in the cancellation of the degeneration of the eigenfunction for a ^ a cr ;t. Consequently, the 
spectrum becomes more complex than in case of no inner core or with uniform conductivity. 



Equatorial symmetry brea king A coupling be t ween d ipole and quadrupo le is the basis of the low di- 
mensional reversal model of Pet relis and Fauvd (|2008l ). iPetrelis et al\ (|2009h where the authors conclude 
that the reversal rate is constrained by breaking of the equatorial symmetry. With this motivation in the 
background, we investigate the influence of equatorial symmetry breaking of the a-effect. In the model an 
additional term is added proportional to cos(2$) leading to 



a = «o sin 



2tt (r - R h 



R 



out 



Rh 



(cosi? + 6cos(2tf)) 



[ID 



1 n characterizes the degree of a spherical Bessel function of a particular free decay mode and determines the radial behavior. 
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where & C 1 and again i? m = 0.35. The main impact of the parameter b is a coupling of dipolar-like and 
quadrupolar-like modes so that a distinction in terms of parity with respect to the equator is no longer 
possible. 

The resulting growth rates for various values of the parameter b are shown in figure [6] whereby a is 
restricted to values around the dynamo threshold. Without symmetry breaking two oscillatory eigenfunc- 




Figure 6. Extraction of growth rates vs. «o f° r dynamo models with equatorial symmetry breaking demonstrating the impact of symmetry 
breaking on the degeneration of the eigcnfunctions. Solid curves denote stationary solutions and dotted curves denote oscillatory solutions. 
b = 0.00,0.01,0.05,0.10. The contour plots on the right hand side show the pattern of the (hemispherical) eigenmodes for case (b) at 
ao = 17 (colour online). 



tions that correspond to an oscillatory dipolar-like mode and an oscillatory quadrupolar-like mode cross 
the dynamo threshold around «o ~ 16.9. The coupling of dipolar-like and quadrupolar-like modes results 
in a number of significant changes in the eigenvalues and geometric structure. Here we restrict the dis- 
cussion to the behavior around the dynamo threshold of the leading eigenmode. First of all, we observe 
a reduction of the critical a with increasing b which is attended by a breakup of the degeneration of the 
leading eigenmodes. However, the dominating eigenmodes retain their oscillatory property but due to the 
coupling of both symmetry classes a distinction into dipolar-like and quadrupolar-like eigenmodes is no 
longer possible. Compared to the undisturbed case, a considerable change in the field structure occurs for 
sufficiently strong a. Around the dynamo threshold (and above) the coupling of both classes results in the 
occurrence of hemispherical dynamos that oscillate in time. Since dipolar-like and quadrupolar-like modes 
contribute roughly the same amount to the new coupled mode (recall that the growth rates of both types 
are roughly equal in the case without any perturbation) the coupled dynamo mode is concentrated within 
one hemisphere. Since dipolar and quadrupolar modes have a very similar structure, their contribution 
cancel out in one hemisphere so that magnetic energy is concentrated in the remaining hemisphere. 
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4 Conclusions 



We have re-confirmed in a simple mean-field dynamo model that oscillatory a 2 -dynamos are possible when 
the radial profile of a exhibits a more complex structure. Our model is based on a particular radial profile of 
the a-effect that behaves oc sin( r). The choice of this p rofile is not arbitrary but is motivated by sim ulations 
of rota ting magnetoconvection ( Giesecke et al. 2005bl ) and quasi-linear computations performed by Sowardl 
(Il974n . On a first glance the occurrence of oscillatory solutions seems quite robust. Around the onset of 
dynamo action, growing oscillatory solutions dominate independently of core conductivity (and/or core 
size) or equatorial symmetry breaking. More important for a disappeara nce of oscillatory solut ions might 
be deviations from the ideal sin-profile such as shifts of the zero crossing ( Giesecke et al\ 2005a ) which has 
not been examined here. 

The eigenfunctions of the a 2 -model exhibit a couple of remarkable properties that characterize the 
geometric structure of the eigenfields. In the oscillatory regime the eigenfunctions are dominated by con- 
tributions with higher radial wavenumber whereas the steady solutions essentially are determined by 
eigenfunctions with sparse radial structure. Hence, it is suggestive to look for an indication of a similar 
contribution of higher radial modes during an reversal in three-dimensional MHD si mulations of the geody- 
namo or in field reconstructions from paleomagnetic observations as e.g. executed bv lLeonhardt and Fabian 
(|2007n . 

For a-distributions with perfect equatorial (anti-)symmetry a oc cos$ the dipolar and quadrupolar 
modes remain separated since no interaction between these modes is possible. The eigenmodes with differ- 
ent equatorial symmetry exhibit an approximate degeneration which would be exact in case of perfectly 
conducting boundary conditions. The corresponding proximity of growth rates for dipole and quadrupole 
has a dramatic impact when a small perturbation is considered that breaks the equatorial symmetry. In the 
vicinity of the dynamo threshold the resulting coupling between dipolar- and quadrupolar-like modes leads 
to hemispherical dynamo action. Indeed, the spatial reconstruction of the last reversal (the Matuyama- 
Brunhes transition ~ 780 kyrs ago) shows a gro wing contribution of the quadrupolar component during 
the actual reversal (ILeonhardt and Fabianl 120071 ) but there are no hints for hemispherical dynamo action 
taken place in the Earth's core. Nevertheless, hemispherical dynamo action might have been the reason 
for th e non-uniform crust magnetization as a result of the ancient martian dynamo (jLandeau and Aubert 
201 lh . 

Our results differ from the achievement of Gallet and Petrelis (2009) who examined a kinematic a 2 
model using a strongly localized a-effect concentrated in two thin shells. In their model it is the coupling 
between dipole and quadrupole induced by equatorial symmetry breaking which leads to an oscillating 
eigenmode, whereas in our model the coupling between different radial modes (induced by the radial 
variation of a) is responsible for the oscillatory behavior. The simplicity of both models, however, prevents 
a robust conclusion which specific behavior indeed is realized in the geodynamo. 

Regarding the Earth's magnetic field, a reason for the suppression of the quadrupole might arise from 
anisotropics of the turbulent flow essentially caused by the fast rotation of the Earth. These anisotropics 
resulting from the fast rotation of the Earth which suppresses variations/fluctuations along the rotation axis 
should be reflected in more realistic models of the a-tensor as well as in the diffusivity tensor. Using more 
realistic anisotropic structure of the a-ten sor (but isotropic rj) results in domination of non-axisymmetric 
(i.e. m = 1) modes ( Giesecke et a/.ll2005af ). This domination can be circumvented by assuming anisotropic 
description for r\ as well (jTilgner 120041 ). Note that non-axisymmetric modes indeed seem to be important in 
the reconstruction of the field pattern during the last reversal (the Matuya ma-Brunhes transition) wher e 
the most important field contribution is determined by a m = 2 contribution ( Leonhardt and Fabian 20071 ). 

A substantial limitation of the presented results is the linear character of the kinematic models, in 
particular regarding the proximity of the eigenvalues of leading dipolar and quadrupolar eigenmodes. 
Thus a robust prediction on the behavior in a saturated state is difficult and more elaborated mean 
field models will require the consideration of backreaction of the magnet ic field by virtue of a - quenc hing. 
Qualitatively, we expect a similar behavior as it has been described by Stefani and Gerbethl ( 20051 ) in a 
one-dimensional non-linear model of a reversing a 2 -dynamo. Assuming that the magnetic quenching of the 
a-effect determines the instantaneous growth rate the authors showed that the corresponding eigenmode is 



February 20, 2012 3:2 Geophysical and Astrophysical Fluid Dynamics giesecke'osc'alpha2 
10 REFERENCES 

inevitably driven towards the oscillatory branch. This shift is attended by self accelerating field decay and 
the emerging local maximum of the growth rate in connection with noise provide the essential preconditions 
to switch from a stationary branch to an oscillatory branch and vice versa. 
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